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Abstract. From random matrix theory it is known that for special values of the 
coupling constant the Calogero-Moser (CM) equation system is nothing but the 
radial part of a generalized harmonic oscillator Schrodinger equation. This allows 
an immediate construction of the solutions by means of a Rodriguez relation. The 
results are easily generalized to arbitrary values of the coupling constant. By this the 
CM equations become nearly trivial. 

As an application an expansion for (e lTr ^ y -') in terms of eigenfunctions of the 
CM equation system is obtained, where X and Y are matrices taken from one of the 
Gaussian ensembles, and the brackets denote an average over the angular variables. 
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1. Introduction 

There is only a very limited number of many-body systems allowing explicit solutions 
P]. This is why the Calogero-Moser (CM) equation system from the very beginning has 
attracted much interest. In its original version the system Hamiltonian reads 



N ^2 i 



n=l " n>m 



(1) 



The parameter g describes the strength of the pair interaction potential decreasing 
with the square of the distance. The eigenfunctions of H, given already by Calogero in 
his original work [2], may be written as 

^nk(x) = z^ 2 ip nk (r)P k (x) , (2) 



where (3=1 + y/1 + Ag, and 

Z J \%n %m\ i T ^ , (*^n %m) ■ (3) 



n>m n>m 



The CM equation system and the ensemble average in the Gaussian ensembles 



2 



x is a short-hand notation for (xi,x 2 , . . . ,xn). The radial eigenf unctions v 9 nfc( r ) 

read 

Vn k (r) = L?> Qr 2 ) e-^ 2 , (4) 

where Ln\x) is a generalized Laguerre polynomial (n = 0, 1, . . .), and 

& = * + i(tf-2) + |tf(tf-l) (5) 

(there is a misprint in Eq. (2.13) of Ref. j2] for b). Pk(x) is a polynomial of degree k, 
totally symmetric in the variables x n . It is obtained as a solution of 

d 2 ^ 1 9 



/~)/y>2 ^^^^ rp />i f~in* 

n n^m 



P k (x) = 0. (6) 



The eigenvalues corresponding to ip n k(x) are given by 

E nk = 1 -{N - 1) + ^N{N -l) + 2n + k. (7) 

Apart from a shift of the energy zero and differing multiplicities of the energy levels 
this is just the spectrum of an harmonic oscillator. 

An equivalent form for the Hamiltonian, which is more convenient in the present 
context, is given by 

2 n=1 \ J 2 ^ x n — x m dx n 

H is obtained from H by introducing centre of mass variables x n = x n — Y^ n x n , and 
applying the transformation 

H = z^ 2 Hz^-H c . m ., (9) 

where H c m , is the centre of mass Hamiltonian. 

From the very beginning there was the strong suspicion that there must be a 
close connection between the CM system and the quantum mechanics of the harmonic 
oscillator, just from the inspection of the spectrum (JZJ). But it lasted for 30 years, until 
the complete equivalence of the two systems was proven by Gurappa, Panigrahi 0. 
For the proof a quite evolved operator technique was applied. There were many futile 
attempts as well to find expressions for the polynomial Pk{x), until recently an explicit 
construction of the solutions of the CM system has been achieved by Ujino, Wadati 0j. 

In this paper it will be demonstrated that there is a much simpler connection 
between the quantum mechanics of the harmonic oscillator and the CM system. In fact 
for the special cases (3 = 1,2,4 this connection will show up to be completely trivial. 

The key ingredient to this result comes from random matrix theory, where the CM 
equations enter via the calculation of averages in the Gaussian ensembles (see e. g. Refs. 
[HIE!)- Since this background probably is not known to all of the experts engaged in 
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the CM equations, a short recapitulation of the aspects relevant in the present context 
is given in Section In the subsequent Section |3] the equivalence of the CM system 
with the Schrodinger equation of the harmonic oscillator is established for the cases 
(3 = 1,2,4. This result is generalized to arbitrary values of (3 in Section HJ In Section 
Elan application to the calculation of averages in the Gaussian ensembles is presented. 
The paper ends with a short discussion in Section |3 



2. The averaging problem in the Gaussian ensembles 

An issue of central importance in random matrix theory is the calculation of averages 
over products and ratios of spectral determinants, 

M v {E i ) = lj\\E i -H\ v \ , (10) 

where the matrix elements of H are assumed to be Gaussian distributed, and the 
Vi may take positive and negative integer values. The standard way to calculate such 
averages uses supersymmetry techniques. In this approach after some steps an integral 
of the type 

I(Y) = [ d[x}^ XY ^™ (11) 



is met, where X, Y are N x N supermatrices. The integral is both over the symmetric 
and the antisymmetric components of X (see Ref. jTj for more details). Diagonalizing 

X, 

X = RXdR- 1 , (12) 
where X D = diag(xi, . . . , x^), one has 

^ (13) 



I{Y) = j (jldxie-^A B(x)f(x,y). 



B(x) is the radial part of the functional determinant, and f(x,y) is given by 

f(x,y) = (eP^) , (14) 

where the brackets denote the average over the angular variables. From the invariance 
property of the trace is is obvious that f(x,y) depends only on the eigenvalues of X 
and Y, but not on the respective angular variables. 

The practical applicability of supersymmetry techniques relies on convenient 
expressions for f(x,y). For Hamiltonians taken from the Gaussian unitary ensemble 
(GUE), the angular average has been calculated by Itzykson and Zuber jH] for ordinary 
matrices. Their result has been generalized to supermatrices by Guhr 0. Explicit 
expressions for the average (jl())l over spectral determinants have been given by Fyodorov, 
Strahov P2|. 
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For the Gaussian orthogonal ensemble (GOE), much more important from the 
practical point of view, and the Gaussian symplectic ensemble (GSE) the situation 
is unsatisfactory. Expressions for the ensemble average of the negative moments of 
the spectral determinant have been given by Fyodorov, Keating jTTj. Brezin and 
Hikami achieved to express averages over products of spectral determinants in terms 
of complicated integrals over determinants of quaternionic matrices Recursive 
relations for f(x,y) with respect to the rank N have been given by Guhr and Kohler 

For all ensembles there are expansions of the type 

fi^^^Pn^Pn^V), (15) 
n 

where the Pn^(x) are Jack polynomials f3J E3- The sum is over all partitions 
n = (m, . . . ,njy), where |n| = J2i n i = ^ is the degree of the polynomial. The 
polynomials are totally symmetric in the variables, and are orthogonal on the interval 
[—§,§] with the weight function ]~f \x n — x m \ , 

n>m 

1/2 1/2 

f dan • • • j dx n ]J \x n - x m f Pl^(x)P^(x) ~ 5 nm , (16) 

-1/2 -1/2 n>m 

see e.g. Ref. jTB]. 

Unfortunately, the results give in the mentioned references are not very convenient 
for the practical application. This is why for the GOE, in contrast to the GUE, only a 
quite small number of analytic results on ensemble averages is available. On the other 
hand there are a number of experiments, in particular in chaotic microwave cavities 
performed in the author's group, where exactly such averages are needed [T71 [TH] - This 
was the original motivation for this work. 

3. The Calogero-Moser equations for {3 — 1,2,4 

f(x,y), considered as a function of x, obeys the differential equation 

(A + TrF 2 ) /(x,?,) = 0, (17) 



where 



and 



TrF 2 = ^y^, = ^^, (18) 



Tl 



is a generalized Laplace operator. Since f(x, y) is independent on the angular variables 
of X, we may restrict the Laplace operator to its radial part. It is well-known from 
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numerous papers on random matrix theory (see e. g. the references cited above), how to 
separate the Laplace operator into its radial and angular part, but for the convenience 
of readers not familiar with the topic the essential steps shall be repeated. 

To keep the discussion simple we assume that all operators are Hermitian, i. e. both 
summation indices i,j in equations (jl8)l and (fTH|) run from 1 to N. The generalization 
to the two other symmetry classes is straightforward though some care has to be taken 
to avoid double counting. For symmetric operators e. g. in all traces the summations 
have to be restricted to i < j. In addition we restrict ourselves to ordinary commuting 
variables. The extension to supersymmetric variables, just as it was demonstrated in 
Ref. (Hj for the GUE, is straightforward as well. 

We start with the somewhat more general equation 

X - (-A + TrX 2 ) i/j n (X) = E n if; n , (20) 

which is nothing but the Schrodinger equation for an independent superposition of 
harmonic oscillators. Its solution may thus by be expressed in terms of products of 
single harmonic oscillator eigenfunctions, or, using Rodriguez' formula for the Hermite 
polynomials, as 



ip nil ...n NN (X) ~ e 



T TrV 2 



n (- v ^) r 



n ( x « ~ v y) 



e 2 



-TrX 2 



TrX z 



(21) 



The corresponding corresponding eigenenergies are 



(22) 



The sub-set of solutions, depending on the radial variables only, is obtained from 
Eq. (|2I|> as 



(23) 



where n is an abbreviations for (ni,ri2, • ■ ■), and 



A* = Tr 



(X - V) k 



rad 



(24) 



is a generalized creation operator. The eigenenergies corresponding to if) n (x) are given 
by 



E n = ^ni + — = \n\ + — 



(25) 
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where \n\ = Y2i n i- Each partition (ni,ri2, ••■) with the same \n\ thus yields an 
independent eigenfunction to the same eigenvalue. 

Only by using invariance properties of the trace we achieved to construct the 
radially symmetric solutions ()23|) of the generalized harmonic oscillator Schrodinger 
equation (f2T)j) . These functions at the same time must be solutions of the radial part 
of the Schrodinger equation. Thus the radial part of the generalized Laplace operator 
entering equation (J2*U|) has to be determined. For the eigenfunctions in addition we need 
expressions for the radial part of Tr(X — V) fc . 

Details of the calculation can be found in | Appendix A For the radial representation 



of the generalized creation operator we obtain 

4 = 0^-^L m > (26) 



where 



D, 



f d | gyv 1 \ (3 1-5 nm ^ 
\ dx n 2 x n Xj~ J 2 x n x m 



Equation (|27J) holds for all Gaussian ensembles with (3 = 1,2,4 for the GOE, the 
GUE, and the GSE, respectively. For the radial part of the Laplace operator we get 

a - = E^^EfV^- (28) 

n n=trn 

Note that this is exactly the operator entering equation © for the polynomials 
Pk(x). For the radial part of the generalized oscillator Schrodinger equation (|2llj) we 
have 



Hip n (x) = E n i, n {x) , (29) 



where 



N 



TT 1 ^ / d 2 2 \ p ^ 1 d 



n=l \ n / n _ Lm n 

This is identical with the CM Hamiltonian (JHJ). 

The derivation of this section has shown that for the special cases f3 = 1, 2, 4 the CM 
equation system is completely trivial. It is nothing but the radial part of the generalized 
harmonic oscillator Schrodinger equation (|20j). As a consequence the solutions are also 
obtained in a trivial way just by picking out all linear combinations of products of 
Hermitian polynomials which are totally symmetric in the variables. 
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4. The general case 

For f3 = 1,2,4 the problem is thus completely solved. But all results can immediately 
be transferred to arbitrary values of (3. 
First we notice that 



^o(x) ~ exp -- V x 2 n ) (31) 




is the ground state eigenfunction of the Calogero-Moser equation (|2~9j) for all values of 
(3, with a ground state energy given by 

*.£ + /J W^i). (32) 

Second, the generalized creation operators obey the commutation rule 

HA\ = A k + {H + k), (33) 

which is the generalization of a well-known relation for the one-dimensional harmonic 
oscillator (see Appendix Bfl . For (3 = 1,2,4 this follows trivially from equation 



But the rule holds for arbitrary values of (3, if only A k + is calculated from expression 
(J2Hj), which is well-defined for all values of /3, and not from equation (J2H), which looses 
its meaning in the general case. 

Since nothing else but this commutation rule is needed for the generalized Rodriguez 
relations to hold, equation (|2~3j) can still be used to generate the eigenfunctions. The 
corresponding spectrum is thus not altered by varying (3, apart from a shift of the ground 
state energy, a fact noticed already by Calogero [2]. 

5. A series expansion for ^e lTr ^ XY ^) 

Generalizing an expansion of e lxy in terms of ordinary oscillator eigenfunctions we are 
now going to expand f(x,y) in terms of the totally symmetric solutions ()23)1 of the 
generalized Schrodinger equation (|2U|) . 

f(x, y) = <e^ y )> = Y, CnmM^Uv) ■ (34) 

nm 

In a sequence of self-explaining steps, including one integration by parts, we have 



e iTr(XY) 



d[Y]ij n (y)f(x,y)= J d[Y]Mv) <e iTr(xy) > 

= c n I d[Y] HTr(Y-\7 Y Pe-^ Y2 

L i 

= c n J d[Y>^ Try2 J] Tr (Y + VyT e iTr(Xy) 
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c n r 



"I 



l[Tr(X-V x P 

i 

l[Tr(X-V x P 



d[Y}e 



-iTrY 2 e iTr(Xy) 



(27T)-il n l^(x) 



(35) 



where c n is the normalization facto of the ip n (x). The differential is given, either in 
Cartesian or radial-angular coordinates, by 



d[Y] = H dY l3 = H\m- y/ J] d v* dn t 



(36) 



i<j 



where dQp is the angular part of the differential. Comparison with equation (}34"|) yields 

c nm = (27r)^iW(p- 1 ) nm , (37) 
where p is the matrix with the elements 

Pnm = J d[X]lp n (x)llj m (x) . (38) 

Eigenf unctions belonging to different \n\ are orthogonal, but for degenerate 
eigenfunctions this is not necessarily the case. The matrix of expansion coefficients 
c nm is thus block-diagonal, where each block corresponds to a given value of \n\. 

The problem of diagonalization of the eigenfunctions belonging to the same \n\ has 
been solved by Ujino, Wadati in a series of papers (THl EI3 EI] ■ The authors introduced a 
new type of totally symmetric polynomials jn (x) they called hidden Jack polynomials 
(or in short, Hi- Jack polynomials), which are obtained from linear combinations of the 
ip n (x) introduced in equation (|2^)) . The Hi- Jack polynomials obey the orthogonality 
relation 



dxi 



d'\\ J { 



x k - x/^exp 



k>l 



k J 



(39) 



Collecting the results we obtain an expansion of f(x, y) in terms of Hi- Jack 
polynomials, 

f(x,y) = <e iT ^)> = (2^ £iH^)^) , (40) 

n 

where the sum is over all partitions n, and 

€^) = c^ 3 ^{x)e^(-W^ \ ■ ( 41 ) 
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c n is a constant normalizing the integral (}39|) to one for n = m. 

Equation (|4"U|) is an alternative to expansion (|15j) in terms of ordinary Jack 
polynomials. It has still to be checked whether this will a real progress in the calculation 
of ensemble averages in the GOE. The formulas given in reference ^H] for the Hi- Jack 
polynomials are quite complicated, but is does not seem hopeless that a more direct 
generation of the polynomials is possible, if the techniques applied in the present work 
are used. 

6. Discussion 

Using the fact that for special values of the coupling constant the CM equation system 
is nothing but the radial part of a generalized harmonic oscillator Schrodinger equation, 
the solution could be immediately constructed in terms of totally symmetric linear 
combinations of products of Hermitian polynomials. The generalization to arbitrary 
values of the coupling constant then was straightforward. 

None of these results is really new. The equivalence of the CM equation system 
with the Schrodinger equation of the harmonic oscillator was suspected from the very 
beginning [2] and was proven recently [B]. Explicit solutions of the equation system, 
too, have been known for some time jlj. In fact the solutions (|23|) derived in the present 
work are completely equivalent to the ones given in Ref. [I]. The fact that the radial 
part of the generalized harmonic oscillator Schrodinger equation reduces to a special 
case of the CM equations, too, has been known from random matrix for several years. 

But the experts working on the CM equations and on random matrix theory did not 
know much of each other as it seems. This is probably why it remained unnoticed for 
such a long time that by a combination of ingredients from both sides the CM equations 
become nearly trivial. 

The results have been applied to derive a new expansion of (e lTr( - xy ^) in terms of 
Hi- Jack polynomials, see equation ([40)1 . Whether this expansion is better adopted to 
random matrix problems as the existing expansion (|T5|) in terms of Jack polynomials 
has still to be explored. 
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Appendix A. The radial part of the momentum operator 

For Hermitian matrices Eq. ()12|) reads, written in components, 




(A.l) 



n 



The CM equation system and the ensemble average in the Gaussian ensembles 10 
Introducing the notation 

J x k , k < N , 

Vk = \ h ^ AT A ' 2 

I ot k , k> N , 

where the a k are the N(N — 1) angular variables, we have 

doc** 



dyk 



where 



and 



, _ i $knO~km ; W j / A /I \ 



Sk,nm — £>k,mn — ^ ] r in ■ (^-5) 



It follows 



hm 

This may be written as 

V,, = ^r m V nm r* m , (A.7) 



where 



9 
<9x„ 



(A.8) 



Using representation (|A.7J) for the matrix elements of V, the radial part of 
Tr(A" — V) fc is now easily calculated. We obtain after a number of elementary steps, 

4 = E[(^- D )*L- ( A - 9 ) 



where 

Equation (jA.lOJl holds for Hermitian matrices, i. e. matrices taken from the unitary 
ensemble. Generalizing the calculation to the other ensembles we obtain equation (|27|). 
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Appendix B. Proof of relation (|33|) 

The D nm and A ra d (see equations (J27)) and (J2HJ) obey the commutation rule 



[^rad, Dnm] = P / J 



\T)nk T) nm ) 2 
\Xk X m ) 



1 



(b.i; 



{x n Xk) 

as can be verified directly from the definitions. The relation is true both for n = m and 
n 7^ m. Introducing the matrix S with elements 



1 &nm 



x \- L 



[Xn X m ) , v"n x k) 

equation (|B.1|) may be written more concisely as 
[A ra d, D nm ] = f3[D, i5] nm . 
Further commutation rules following directly from the definitions are 



(B.2) 



(B-3) 



[A ra d, Xn] 

x ki D T , 



-2D r 



^^nmXn 



Introducing the operator A + with the components 



Xn^nm D nm , 



(A- 



one obtains by combining equations (jB.3J) to (jB.5J) 

(A+) n J = ( A +)nm ~ 7} [ A +i S ]nm ■ 

Now the wanted commutator ([33)1 can be calculated, 



(B.4) 
(B.5) 

(B.6) 
(B.7) 



[H,A k + ] =J2\ H ^ A 



l +J J urn 



nm 
fc-1 



i/=0 ijnm 
P 



v (fc-«/-l)l 



fe-1 



v(fc-«/-l)l 



i>=0 ijnm 



(B-8) 



The second term on the right hand side does not contribute as is easily seen, and 
we are left with equation q. e.d. 
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